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ABSTRACT 

Context. Magnetic fields are believed to play a crucial role in the process of star formation. 

Aims. We compare high-angular resolution observations of the submillimeter polarized emission of NGC 1333 IRAS 4A, 
tracing the magnetic field around a low-mass protostar, with models of the collapse of magnetized molecular cloud 
cores. 

Methods. Assuming a uniform dust alignment efficiency, we computed the Stokes parameters and synthetic polarization 
maps from the model density and magnetic field distribution by integrations along the line-of-sight and convolution 
with the interferometric response. 

Results. The synthetic maps are in good agreement with the data. The best-fitting models were obtained for a protostellar 
mass of 0.8 Mq, of age 9 x 10* yr, formed in a cloud with an initial mass-to- flux ratio ~ 2 times the critical value. 
Conclusions. The magnetic field morphology in NGC 1333 IRAS 4A is consistent with the standard theoretical scenario 
for the formation of solar-type stars, where well-ordered, large-scale, rather than turbulent, magnetic fields control the 
evolution and collapse of the molecular cloud cores from which stars form. 
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■ 1. Introduction 



NGC 1333 IRAS 4A (herafter IRAS 4 A) is one of the 
prototypical low-mass young stellar systems in the earli- 
est stages of evolution, and it is still deeply embedded in 
an infalling dense molecular and dusty envelope (Sandell 
et al. 1991; Di Francesco et al. 2001) and powering a well- 
colHmatcd outflow (Blake et al. 1995; Choi 2005; Choi et 
al. 2006). The BIMA spectropolarimetric observations have 
detected and partially resolved the polarization in both the 
dust (at 1.3mm) and line (CO J = 2-1) emission (Girart 
et al. 1999). Recent observations with the submillimeter 
Array (SMA) at 877 /im (see Fig. [Ij show that the mag- 
netic field associated with the infalling envelope has a well- 
defined hourglass morphology on scales of a few hundred 
AU (Girart et al. 2006, hereafter GRM06). In this paper, 
we perform a quantitive comparison of the observed sub- 
millimeter polarization data with models of the collapse of 
magnetized molecular cloud cores and show that the data 
support the theoretical scenario where the ordered, mean 
component of the interstellar magnetic field controls the 
evolution and collapse of the molecular cloud cores from 
which stars form (see e.g. Shu et al. 1987, 1999; Machida 
et al. 2005a,b; Banerjee & Pudritz 2006). 
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Fig. 1. Map of NGC 1333 IRAS 4A, from GRM06. 
Contours show the continuum emission at 877 /im, bars 
indicate the direction and degree of polarization (magnetic 
field vectors), and the color map shows the polarized inten- 
sity. At the distance of 300 pc, 1" corresponds to 300 AU. 
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2. The models 

We have adopted two models for the magnetic field 
in of IRAS 4A. The first is an axisymmetric, time- 
dependent, calculation of the inside-out collapse of an uni- 
formly magnetized cloud with ambipolar diffusion (Galli 
& Shu 1993a, b, hereafter GS93). The second is a steady- 
state, axisymmetric calculation of the accretion flow on a 
low-mass protostar and the associated magnetic field, in- 
cluding the effects of Ohmic dissipation (Shu et al. 2006, 
hereafter S06). Both studies ignore the effects of rotation 
on the collapse (for this, see Galli et al. 2006). 

The study of GS93 focused on the formation of a 
large flattened structure ( "pseudo-disk" ) around an accret- 
ing protostar due to the effect of the strong non-radial 
component of the Lorentz force on the collapsing matter. 
According to GS93, ambipolar diffusion during the dynam- 
ical collapse phase plays a minor role on scales larger than 
a few 10^ AU, especially at very early times after the on- 
set of collapse. Under quasi-fleld freezing, the infalling gas 
pulls the magnetic fleld towards the center, strongly pinch- 
ing the fleld lines near the protostar. The parameters of 
the GS93 model are the spatial scale of the core prior to 
collapse, rp , and the nondimensional time elapsed since the 
onset of collapse, r. The scale ro is defined in terms of the 
effective sound speed Ues (including thermal, turbulent and 
magnetic support) and the initial magnetic field Bq by the 
expression rg = 2a^ff/(G^/^i?o), where G is the universal 
gravity constant. The physical conditions in the core that 
gave birth to IRAS 4A are difficult to derive. In particular, 
no Zeeman measurement of the magnetic field strength in 
the surroundings of IRAS 4A is available. Thus, we adopted 
the "fiducial" values of GS93, a^s — 0.35 km s""'^, and Bq — 
30 /iG. Then ro = 2.1 x 10'* AU, and the nondimensional 
time r is given by r = t/tQ with Iq = ro/ucs = 2.9 x 10^ yr. 
The results can be easily scaled to different values of tq. 

The only parameter of the S06 model is the Ohm's ra- 
dius, a spatial scale associated with Ohmic dissipation and 
defined as rohm — '7^/(2G'M^,), where rj is the Ohmic resis- 
tivity, assumed to be spatially constant. With a non-zero 
resistivity, the magnetic field lines strongly bent inward by 
the pull of the infalling gas relax to an almost straight and 
uniform field configuration in a region of a size ~ rohm, with 
the enclosed magnetic flux reduced with respect to the field- 
freezing value (by a factor of ^ 100 at r = rohm)- The value 
of T] depends on the ionization, temperature, and chemical 
composition of the gas close to the source and is not well- 
constrained. Moreover, the value of 77 is probably associated 
to macroscopic plasma instabilities rather than microscopic 
collisional processes (Shu et al. 2007). Estimates of rohm are 
in the range 1-100 AU for a solar- mass star (SOB). 

3. The method 

In spite of the idealization of the models, a comparison 
with the data is not straightforward. First, the orientation 
in space of the models is specified by two viewing angles, the 
position angle of the polar axis with respect to a reference 
direction in the plane of the sky (with = corresponding 
to a polar axis aligned north-south), and the inclination ip 
of the meridional plane of the model with respect to the 
plane of the sky (with = corresponding to a polar 
axis in the plane of the sky). These orientation parameters 
are well-constrained by the observations. In addition, the 



GS93 model depends on two parameters, a spatial scale 
ro and time, whereas the steady-state S06 model depends 
on a single parameter, the Ohm's radius rohm- The latter 
in particular is not well-constrained by the theory (S06). 
Thus, the data must be compared with an extensive set of 
model density and magnetic field distributions obtained for 
different values of the model parameters and the orientation 
angles (j) and ip. 

First from a grid of density and magnetic field models 
we obtained a large number of maps of Stokes parameters /, 
Q, and U, using the same method as described in Gongalves 
et al. (2005). Once the models were generated, they were 
convolved with the SMA interferometric response. This was 
done by converting the modeled map to visibilities using the 
same distribution of visibilities in the u, v plane and the 
same u, v weighting as the SMA observations of IRAS 4A. 
This procedure was followed independently for the Stokes 
/, Q, and U. Once the synthetic SMA-like Stokes /, Q, U 
maps were obtained, the polarization intensity and position 
angle were obtained in the same way for the SMA data 
and for the models. Since this convolution is nonlinear, it 
is necessary to scale the model intensity with that of the 
observations, which is equivalent to selecting a specific value 
for each map for the product kB{T), or, in other words, a 
mass calibration factor. This is performed with an iterative 
procedure. 

The third and last step was to compare the modeled 
and observed polarization angle at each observed position. 
A quantitive evaluation of the goodness of the fit is given 
by the standard deviation of the distribution of the resid- 
uals between the modeled and the measured polarization 
angles. The continuum emission and the polarized inten- 
sity predicted by the model were also used as additional 
constraints. 

4. Results 

4.1. The G593 model 

We have computed a total of 270 maps for r = 0.3, 0.5, and 
0.7, with different values of (p and -0- The best-fit models 
have (j) = 50°-60°, ip = 0°-30°, ro = 3.2 x lO" cm, and 
T = 0.3. Although the fit is not very sensitive to the value of 
f/), the best results are obtained for small inclinations of the 
model with respect to the plane of the sky, with the pseudo- 
disk seen almost edge-on. In Fig. [2] we show the histograms 
of the residuals in polarization angles (differences between 
predicted and observed values) for the three values of r. 
The standard deviations are a — 14.8°, 15.6°, and 18.5° for 
T — 0.3, 0.5, and 0.7 (the instrumental uncertainty on the 
polarization angles is Uinstr. = 6.2°, see GRM06). 

Although the polarization angles and the intensity pro- 
file of IRAS 4A can also be satisfactorily reproduced by 
models with r — 0.5-0.7, the isocontours of the emitted flux 
are flatter in these cases than observed, and the polarized 
flux intensity not reproduced as well as in the r = 0.3 case. 
For T = 0.7, neither the shape of the isodensity contours 
nor the polarized intensity distribution can be reasonably 
reproduced for any possible viewing angle. In Fig. [3] we 
show a comparison of the modeled and observed intensity 
and polarization map for the best-flt value r = 0.3. For this 
case, the largest contribution to a comes from a systematic 
deviation of the polarization vectors on the top-left side of 
the map. 
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Fig. 2. Histogram of residuals in polarization angles for 
the GS93 models. Solid line, r = 0.3; short-dashed line, 
r = 0.5 AU; long-dashed line, r = 0.7 (histograms slightly 
displaced for clarity) . The dotted curve shows the Gaussian 
distribution of the measurement uncertainty with CTinstr = 
6.2° (GRM06). 



With our choice of parameters, the best-fit model with 
T — 0.3 implies an age for this protostar oft* = 8.6 x 10'' yr 
and a mass of « ajJfjt*/G = 0.8 Mq. The mass of the 
infalling envelope within a radius r = 10'^ AU from the 
source is 0.5 Mq, in agreement with the mass distribu- 
tion derived by Belloche et al. (2006) from interferometric 
and single-dish continuum observations (0.7 and 0.9 Mq, 
respectively). From the best- fit model, the magnetic flux 
enclosed in a radius r ~ 10^ AU is $ = 2.4 x 10^° G cm^, 
so the mass-to- flux ratio of the envelope plus star is ^ 1.7 
times the critical value (M/$)cr = {2ttG^^^)~^ , in agree- 
ment with the independent estimate of GRM06. Since the 
effects of ambipolar diffusion are negligible on the collapse 
time scale in the GS93 model, this value characterizes the 
mass-to- flux ratio of the original 1.3 Mq of cloud core that 
have collapsed to form the protostar and the envelope. 



\ 


\ 


\ 


1 


1 


\ 


1 








\ 


\ 


\ 


\ 


\ 




^ 1 






\ 


\ 


V 




'T 


\ 




\t 






\ \ 






\ 


\ 


\ 




A 


\ \. 


K 


\ 


\ 






J 


/ 




1 


\ \( 


\ 


\ 






\ 




V; 




- 






















\ 1 




















\,/ 




















■— . 








K 












— \ 


— 
































/\ 






/ 


1/ 








V 








\ " 


1 




/ 








\ \ 




N \" 


1 




J 


1 


\ 


\ 


\ V 




\ \ " 




1 






\ 






^\ 


\ 


\ 






1 


1 


\ 


\ 


\ 


\ 


\ 







arcsec 



-2 



Fig. 3. Synthetic maps for the GS93a model: continuum 
emission (contours), polarized emission (greyscale), polar- 
ization degree (segments) . Red segments show the observed 
polarization vectors. Model GS93 with r = 0.3, ^ = 30°, 
6 = 55°. 



tern alone is not sufficient for constraining the model pa- 
rameters. Models with high values of rohm (''ohm ~ 200- 
500 AU) are not able to reproduce the observed polarization 
pattern very well or even the two-lobe structure of the po- 
larized continuum emission. The best fit to the polarization 
angles is obtained for models with rohm ^ 50 AU, suggest- 
ing that the actual value of rohm may be close to or below 
the spatial resolution of the observations (^ 10^ AU). In 
fact, the two-lobe shape of the polarized intensity emission 
evident in Fig. [T] can only be reproduced with very low 
values of the Ohm's radius, rohm = 5-50 AU, suggesting, 
in agreement with the findings of Sect. 14.11 that the mag- 
netic field down to the scale of resolution is dominated by a 
pinched, rather than uniform, component. Figure |4] shows 
the histograms of the polarization angle residuals for the 
models with rohm = 5, 50 and 500 AU, characterized by 
standard deviations a = 11.7°, 13.8°, and 22.6°. The po- 
larization and intensity map for the case rohm = 5 AU is 
shown in Fig. [5l 



4.2. The S06 model 

The S06 model is only applicable to a small region with re- 
spect to acfft (~ 6000 AU with the parameters of Sect. 14. 1| ). 
and assumes steady-state over this region. Thus, no ages or 
stellar masses can be derived for this model. In principle, a 
comparison of the polarization data for IRAS 4A with the 
S06 model is able to constrain the value of rohm without 
relying on uncertain estimates for the Ohmic resistivity. 

For the S06 model, we computed a total of ~ 450 syn- 
thetic maps with different values of (j) and and values of 
'''Ohm ranging from 5 to 500 AU. Best-fitting models are 
obtained with = 60°-65° and V = 25°-55°. The results 
for the GS93 model discussed in Sect. 14. II are not very sen- 
sitive to the value of "0, and the observed polarization pat- 



5. Discussion and conclusions 

From the discussion of Sect. [H it is clear that the agree- 
ment of magnetic collapse models with the continuum po- 
larization data for IRAS 4A is very good, supporting the 
standard theoretical scenario for the formation of low-mass 
stars from magnetized molecular cloud cores. The residuals 
distributions indicate that the best agreement with the data 
is obtained with GS93 models with early ages (ct — 14.8° 
for T — 0.3) or for S06 models with low values of the Ohm's 
radius (cr = 11.7° for rohm — 5 AU). The low value of 
''Ohm obtained from the best-fitting S06 model (rohm = 5- 
50 AU) might also indicate that the steady-state conditions 
assumed in that model have not yet been reached, and the 
size of the region of strong magnetic dissipation is below the 
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Fig. 4. Histogram of residuals in polarization angles for 
the S06 models. Solid line, rohm = 5 AU; short-dashed 
line, rohm = 50 AU; long-dashed line, rohm = 500 AU 
(histograms slightly displaced for clarity). The dotted curve 
shows the measurement uncertainty as in Fig. [2] 
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Fig. 5. Synthetic maps for the S06 model: continuum 
emission (contours), polarized intensity (contours and 
greyscale), polarization degree (segments). Red segments 
show the observed polarization vectors. Model with rohm = 
5 AU,^ = 54°, (j) = 65°. 



resolution limit of the observation. At any rate, both mod- 
els suggest that the magnetic field lines in a region of radius 
~ 500 AU from the central source(s) are almost radial, as 
expected for collapse under ideal MHD conditions. 



The almost radial geometry of magnetic field lines in 
IRAS 4A is also supported by the distribution of polarized 
intensity shown in Fig. [1] In fact, only models with strong 
central concentration of magnetic field lines reproduce the 
observed two-lobe structure of the observed polarized in- 
tensity, with two large "polarization holes" in the midplane 
of the pseudo-disk, on both sides of the central protostar. 
Rather than a variable efhciency alignement of dust grains, 
the presence and extent of these "holes" seem to indicate 
a combination of geometrical effect (many field lines point 
to the observer in a strongly pinched geometry) combined 
with beam dilution effects (the pinched region is smaller 
than the beam), resulting in a cancellation of Stokes pa- 
rameters along the line-of-sight. 

Some amount of field dissipation must have already oc- 
curred in the central regions of IRAS 4A. The presence 
of two continuum peaks in IRAS 4A separated by 1.8" 
(Looney et al 2000; Reipurth et al 2002), corresponding 
to 540 AU at the assumed distance of 300 pc, suggests 
that some fragmentation of the pseudo-disk has already 
occurred, leading to the formation of a (possibly bound) 
protobinary system. Without field dissipation, the catas- 
trophic magnetic braking associated to the concentration 
of magnetic fields in ideal-MHD collapse calculation pro- 
vides a fierce opposition to fragmentation and to the for- 
mation of circumstellar disks (Galli et al. 2006; Hennebelle 
& Fromang 2008; Mellon & Li 2008). 

From the distribution of residuals of the fit, it is pos- 
sible to estimate the relative magnitude of the turbulent 
vs. the ordered component of the magnetic field using the 
formula \SB\/\B\ « (tr^ _ af^^^^y/^ (GRM06). With the 
a obtained for the best-fitting GS93 and S06 models, we 
derive |(5B|/|B| « 20%. This is probably an upper limit on 
the intensity of the turbulent field. Thus, The magnetic field 
morphology in IRAS 4A is consistent with the standard the- 
oretical scenario for the formation of low-mass stars from 
cores threaded by ordered rather than turbulent magnetic 
fields. 
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